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Numerical  Analysis  of  a  Singular  Integral  Equation 
Arising  from  Electromagnetic  Interior  Scattering 


1  Introduction 

For  safety  and  health  reasons,  it  is  of  considerable  interest  to  assess  the  short-  and  long-term 
effects  of  electromagnetic  (EM)  radiation  on  people  working  near  radars  and  other  similar 
EM- wave-generating  devices.  Research  to  understand  this  can  be  classified  as  epidemio¬ 
logical,  experimental,  and  numerical.  In  numerical  electromagnetic  dosimetry  one  is  led 
naturally  to  the  problem  of  solving  the  Maxwell’s  Equations  inside  a  highly  inhomogeneous 
and  highly  dispersive  body.  One  of  the  solution  approaches  is  to  solve  an  equivalent  problem 
in  the  frequency  domain  using  a  volume  integral  equation  formulation. 

Mathematically,  in  the  time-harmonic  case,  if  the  body  ( V  )  is  incident  by  an  electric  field 
E*(r)  and  if  E(r)  is  the  total  electric  field  inside  the  body  (r  €  F),  then  the  scattered  field, 
defined  as  the  difference  between  the  two  fields,  Es(r)  :=  E(r)  -E'(r),  can  be  shown  to  take 


the  form 

Es(r)  =  (1+  1  VV-)  jf  g( r,r')  F£(r')  dV' 

(1) 

in  which 

F£(r)  :=  r(r)E(r) 

'  r(r)  :=  k\ r)  -  fc; \ 

pjk0r 

9<r'r)  ;=  4*r 

r  :=  1  r  —  r'  1 

Here  j  —  \[—  1  and  kD  and  k( r)  are  the  wave  numbers  associated  with  free  space  and  the 
body  respectively.  Further  manipulation  of  Equation  (1)  to  move  the  differentiations  under 
the  integral  signs  results  readily  in  a  vector  integral  equation  of  the  form  (see  [7]  for  details) 

Es(r)  =  A(r)FE(r)  +  /  G(r,r')(FE(r')  -  F£(r))  dV'  (2) 

JV 

in  which  the  dyad  A(r)  becomes  unbounded  near  the  boundary  and  the  dyad  (the  dyadic 
Green’s  function)  G(r,r')  has  a  singularity  of  type  0(r~3).  Previous  attempts  to  solve 
this  equation  [6,  7]  using  a  Moment  Method  or  a  Nystrom  Method  have  been  successful  only 
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for  a  restricted  class  of  parameters.  Before  attempting  to  conduct  a  thorough  numerical 
analysis  of  Equation  (2),  we  opted  instead  to  analyze  the  following  simpler  but  analogous 
1-D  integral  equation: 

A  {t)4W-Kg<f>(t)-Kl4>{t)  =  xM  (3) 

for  t  £  [a,  b]  and  x,  £  C([a,b]).  Here,  analogous  to  Equation  (2), 

KM1)  ='  [  K(t,s)<Ks)ds 

**  a 

b 

KJ(t)  =  [  kt{t,s)[<j>{s)-<j>{t)}ds 

J  a 

The  functions  A (t),kb(t,s),  and  kg{t,s)  also  have  properties  analogous  to  that  m  the  3-D 
case7  namely: 

•  X(t)  >  0  on  (a,  b ) 

•  lim  A(£)  =  oo  and  limA(t)  =  oo. 

•  h(t>s)  =  7,  >  0. 

•  kg(t,s)  is  continuous  on  [a,  b\  X  [a,fr]- 

It  should  be  mentioned  that  the  problem  considered  here  is  not  equivalent  to  the  1-D 
Maxwell’s  Equations  wherein  E  is  dependent  on  only  one  spatial  dimension.  It  is  well-known 
that  the  Green’s  function  for  the  1-D  Maxwell’s  Equations  is  much  better  behaved.  To  keep 
the  problem  here  simpler,  however,  we  will  ignore  the  ‘good  kernel  kg(t,s)  in  the  following 
analysis  and  concentrate  only  on  the  equation 

A  (t)4>{t)-Ki<Kt)  =  xAt)  (4) 

As  a  preliminary  analysis  of  Equation  (3),  we  studied  numerically  the  1-D  problem  anal¬ 
ogous  to  Equation  (2)  in  which 

g(t,s)  =  \t  -  s|(ln  \t  -  s\  —  1). 

After  properly  moving  the  double  differentiations  under  the  integral  sign,  the  resulting 
k  kb(t,s )  and  A(t)  can  be  shown  to  satisfy  the  properties  mentioned  above.  Numerical 

solutions  were  obtained  using  several  variants  of  the  Nystrom  methods. 

1.  Product  method  with  extrapolation  at  the  end  intervals 

2.  Gauss-Legendre  method 
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3.  A  simple  method  (in  which  uniformly  spaced  integration  points  and  uniform  weights 
are  used) 

In  each  case,  apparent  convergence  was  obtained.  The  main  result  of  this  study  is  a  rigorous 
mathematical  proof  of  the  convergence  of  a  numerical  method  used  to  solve  Equation  (3). 

In  Section  2,  we  will  re-formulate  the  problem  and  put  it  into  perspective.  In  Section 
3,  we  will  investigate  the  properties  of  an  operator  K  that  arises  from  the  re-formulated 
problem.  The  numerical  method  used  to  solve  the  problem  will  be  defined  in  Section  4, 
and  some  preliminary  properties  of  the  associated  numerical  integral  operators  Kn  will  be 
explored.  Several  convergence  theorems  for  the  numerical  integral  operators  Kn  will  be 
proven  in  Section  5.  In  Section  6,  a  convergence  theorem  for  the  numerical  solution  of  the 
complete  problem  will  be  given.  Finally,  we  will  conclude  with  some  closing  remarks  in 
Section  7. 

2  Statement  of  the  Problem 


We  begin  with  the  integral  equation  in  (4),  namely 

X(t)  (f)(t)  -  ^  f  |f-5|_1[^(s) -<£(«)]  ds  =  Xi(i)  (5) 

*'a 

defined  on  an  interval  (a,  b).  Here  we  assume  the  constant  >  0  and  the  function  Xi  £ 
C([a7b]).  Furthermore,  we  assume  the  function  A  is  positive  and  continuous  on  (a,  b)  and 
that  X(t)  — y  oo  as  t  — y  a  and  as  t  — >  6.  Under  these  assumptions,  Equation  (5)  can  be 
transformed  to  an  integral  equation  of  the  second  kind: 

H1)  -  7(0  /  .  |<  -  -HM  -  <£(01  ds  =  x(t)  (6) 

J  a 

defined  on  the  interval  [a,  b].  Here  7 (t)  :=  7, /A (£)  >  0  on  (a,  b)  and  together  with  x(£)  := 
X,(t)/A(t)  may  be  assumed  to  be  continuous  on  [a,  b] ,  because  of  the  assumptions  on  A. 
Define  the  operator 

K<j>(t)  :=  f  |£  —  •s|~1[^>(s)  —  4>{t)]  ds  (7) 

J  a 

Then  Equation  (6)  can  be  written  in  the  familiar  operator  notation  as 

(I-1K)<f>  =  X  (8) 

Our  problem  is  then  to  analyze  the  numerical  solution  of  this  equation  when,  in  particular, 
a  simple  Nystrom  method  (to  be  described  in  Section  4  below)  is  used. 
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Tlie  problem  being  addressed  here  differs  from  conventional  weakly  singular  integral 
equations  in  at  least  two  fundamental  ways.  First,  while  Equation  (6)  contains  the  difference 
term  used  in  the  well-known  Singularity  Subtraction  Method,  namely 

—  f  k(t,s)[<j)(s)  -  (f>{t)]ds, 

"  a 

the  subtracted  term  b 

J  k(t,s)<j>{t)ds  =  jf  k{t,s)ds 

in  our  case  is  divergent.  This  is  in  stark  contrast  to  the  conventional  case  where  the  sub- 
tracted  term  is  and  must  be  finite. 

Second,  for  weakly  singular  integral  of  the  second  kind 

I  ds  =  x(0 

J  a 

or 

( I-Kw)4>  =  x 

where  |MM)|  <  C  \s  -  t|“_\0  <  a  <  1,  I<w  is  compact  from  C([a,b])  -*  C([a,b]).  Conse¬ 
quently,  the  analysis  of  a  typical  numerical  method  taking  the  form 

(f  -  Kn)<f>„  =  X  (9) 

can  be  based  on  Anselone’s  Collectively  Compact  Operators,  wherein  the  operators  Kn  are 
each  compact  from  C([a,b])  ->  C([a,b]).  (See,  for  example,  [1]).  Unfortunately,  m  our 
problem  the  operators  are  not  compact,  as  we  shall  see  below. 

3  Mapping  Properties  of  K 

We  first  investigate  the  mapping  properties  of  the  operator  in  Equation  (7). 

:=  [  \t  -  s\~l  [^(s)  -  (j>{t)]  ds 

J  a 

We  recall  a  function  f  is  uniformly  Holder  continuous  of  order  a  (0  <  a  <  1)  on  an 
interval  [n,  h]  if  there  exists  a  constant  C  such  that 

I  fix)  -  f{y)  l<  G  I  x~y\a 
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for  all  x  and  y  in  [a,  6].  Define 


{The  space  of  all  uniformly  Holder 
continuous  functions  of  order  a 
on  an  interval  [a,  fc] 

The  following  properties  of  C(0’Q)([a,6])  are  well-known: 

Proposition  2.  For  0  <  a  <  (3  <  1 

1. 

2. 

3. 


4. 

Proof.  (See  [5,  3])  U 

Corollary  3.  For  0  <  a  <  (3  <  1,  the  mapping  from  (<7(  ’^([a, h]),  ||  •  ||J  onto 
(C^Qa.t]),  ||  -  11^)  is  unbounded. 

Proof.  Otherwise,  the  identity  map  1^  would  be  compact  on  the  infinite  dimension  space 

cP*MY  u 

Corollary  4.  For  0  <  oc  <  /?  ^  1,  ( [<x ,  fc] )  is  not  a  Banach  subspace  of  C  ([a,  6]). 

Proof.  Else  \~l  would  be  bounded,  by  the  closed  graph  theorem  [4],  since  a  and  therefore 

I"1  are  closed  operators.  LI 

Proposition  5.  For  0  <  a  <  <  1,  K  :  C<"',([a, 6])  ->  C1"'"'  ([a, 4])  is  compact. 


Cm(M))  C  Cl,-’([«,6]) 

([a,  b])  is  a  subalgebra  of  C([a}  6]) 

6])  is  a  Banach  space  under  the  norm 

ll/IU  =  ll/IU  +  i/I. 
I/M -/fall , 


where 


l/L  :==  sup{- 


\x  —  y|c 


±  y } 


(10) 


is  a  semi-norm. 


Imbedding  maps  1  \  ’^([a,  b])  — >  CK  are  coniPacl 
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proof.  Mimicking  the  steps  in  one  of  the  proofs  in  [3],  one  can  show  that  K  is  bounded  from 
((aj  b\)  ->  C(M)(K  *>]),  where  S  :=  (a  +  j3)/2.  Using  the  fact  that  the  imbedding  from 
C{M)  ([«,&])  — >  C(°’a)  ([a,b])  is  compact,  the  proposition  follows  immediately.  U 


Unfortunately,  classical  Fredholm  theory  does  not  apply  here,  because  of  the  following  ob- 
servation. 

Proposition^.  C'(0'“)  ([a,  b ])  is  not  invariant  under  K  for  any  0  <  a^<  1. 

proof.  By  direct  calculation,  one  can  show  that  <j>{t)  :=  &  €  C(  ‘P>([a,5]),  but  K<f>{t)  <£ 

(7(o’0)([a,  6]),  assuming,  without  loss  of  generality,  [a,b]  =  [0,1]-  LJ 


For  theoretical  as  well  as  numerical  reasons,  it  is  desirable  to  consider  operators  L  whose 
range  is  contained  in  its  domain,  so  that  L2,  for  example,  is  defined.  This  leads  us  to  he 
following  spaces.  For 


0  <  a  <  1,  we  define 


X„:=U{C,"’(M])l°<^  1} 

In  particular.  X,  is  the  set  of  all  functions  defined  on  [a,  b]  which  are  uniformly  Holder 
continuous  of  some  order  a  G  (0, 1]. 


Lemma  7.  For  0  <  a  <  1,  (Xa,  ||  •  [|J  is  a  normed  hnear  space  and  is  invariant  under  K. 
proof,  (X„,||  -  |U)  is  a  linear  subspace  of  C'0'"' ([a,  b}),  ||  -  ][„)-  The  invariance  follows  from 

Proposition  4.  LI 

While  the  semi-norm  |  •  |.  in  Equation  (10)  and  hence  the  norm  ||  •  ||„  are  defined  for  a  6  (0. 1] 
on  £7(0,a) ([a,  6]),  it  is  convenient  (and  also  consistent)  to  define 

ll/llo  :=  ll/IL.  /  G  X0 

Lemma  8.  (X„,  ||  •  ||.)  is  a  normed  linear  space  and  is  invariant  under  K. 

II  -II.)  =  (X.,  II  ■  ID  is  a  linear  subspace  of  (C(M),  II '  ID.  and  the  rnvariance 

follows  again  from  Proposition  4.  U 


6 


While  not  germane  to  our  discussion  here,  it  can  be  shown  that  the  closure ol  t*.,  II 
C(°-a\[a,b\)  is  not  C{°'o){[a,b}),  even  though  &  “’(M)  contains  C  ((a,  b})  for  all  (3  > 


the  closure  of  (X„,  i|  - 1|„)  in 
,  n(0'p)t\n  h}\  for  all  B  >  a. 


Proposition  9.  K  is  unbounded  on  (X„,  ||  •  ||J  f°r  any  0  <  a  <  1.  ^  ^ 

PropE"- Assume,  without  loss  of  generahty,  [a,b]  =  10,1].  One  can  then  readily  show  that 

:=  ^+1/” 

is  a  bounded  sequence  in  Xa,  but  \\K<f>„\\a  ->  oo  as  n  ->  oo.  U 

4  Numerical  Integral  Operators,  Kn 

The  numerical  solution  of  Equation  (8)  can  be  defined  in  terms  of  the  following  nu“er^ 
•  t  i  U  is  basically  the  operators  associated  with  the  Nystrom  method. 

>  0.  1«  on  the  interval  M 

“to  n  e^  subintervals.  We  associate  with  the  partition  P.  the  operator  if.  defined  on 

C([a,  b})  as  follows. 


where 


KM)  ••=  E  w-j  &*»■(*) 

3= 1 


k  =  2n 


(f>  G  C([a,b]) 


wnJ  =  ( b—a)jkn — : 

K 

II 

-si 

t  G  [t„j- 1  >  ^n.,] 

1  VK 

a  „,m  =  wj-m 

t  .  —  a  +  jh  ,  j  —  0,...,kn 

n'}  1 

t*  =  (*»„•_ i+4nj)/2>  i  =  h--A 


For  simplicity  we  have  purposely  chosen  each  weight  wn  .  associated  with  the  j-th  subin- 
terll  P  to  W  dependent  o^ly  on  the  integer  n  and  not  on  j.  More  sophrst.ca^d  cWe 
for  the  weights  is  of  course  possible,  but  the  resulting  analyses  would  be  more  cornea  . 
Also  K  4  is  actually  defined  for  any  function  4  that  is  merely  defined  on  the  nt  (  ,  f 
However  here  we  are  only  interested  in  those  functions  that  are  at  leas  contrnuous.  It 
obvious  that  if  4  €  C([a,b]),  then  so  is  KJ.  That  is,  C([a,b])  rs  rnvanant  under  K„. 
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The  numerical  solution  <f>n  to  Equation  (8)  is  now  obtained  by  solving  the  equation 

(l-iKn)<f>„  =  x  (11) 

By  collocation  at  the  kn  mid  points  {£*  },  the  following  system  of  linear  equations  is  obtained 

[I -lKn)<j>n{t*n  .)  =  x(Cj.)»  3  =  1,  -A  (12) 

from  which  {0rt(£*  ^)}  can  be  solved. 

For  comparison  with  the  operator  K 1  we  will  look  at  some  properties  of  Kn.  Unless 
specified  otherwise,  we  will  always  assume  n  is  a  positive  integer  in  the  following.  Also  for 
later  convenience  we  introduce  the  following  functions,  each  of  which  depends  only  on  n: 

V’nW  =  =  E  W~j9nj(t) 

■J=l 

We  will  first  look  at  some  properties  of  Kn  on  C([a,b])  and  then  its  properties  on 
C{0>ex)  ([a,  6]),  ol  <E  (0,1].  As  we  have  already  noted,  we  have 

Proposition  10.  C([a,b])  is  invariant  under  Kn. 

Moreover,  we  have 


1 

Proposition  11.  Kn  is  bounded  on  C([a,b])  with  =  8  £  jVy 

Proof.  From  the  definition  of  Kny  it  follows  immediately  that 

KJ(t)  =  +  KaW) 


where 

i=i 

K.A‘)  = 

Now  Kn ,  is  compact  and  hence  bounded  on  C([a,b]),  because  it  has  finite  dimensional  range. 
Since  gni(t)  €E  C([a,b]),  so  does  ij>n(t).  Hence  Kn2  is  also  bounded  on  C([a,b]).  It  follows  that 
K„  must  be  bounded  on  C([a,b]).  To  find  the  norm  of  Kn  on  C([a,b]),  let  t*  =  (a  +  b)/ 2. 
One  can  verify  directly  that 

W.IL  =*.(*') 


t 
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For  any  (j)  G  C([a,b]), 


\Kn<f>(t)\  <  2||^|L^(i) 

<  2||^|L^(f) 

Hence  <  2*.(f).  Since  f  4  K,,  1  =  h-  -,K,  ^ere  exists  *. .  6 ^C([»,b])  such 

that  ||*.||_  =  1  and  *.(f)  =  =  1,  )  =  l.-.*.-  Th"s.  * 

H-KJloo  —  2V’n(t*)-  Ifc  follows  that 

IKII.=2iUf)' 

&n_ j  ^n—1  ^ 

By  direct  verification,  one  obtains  =  4  £  spi-  Hence,-||Kn||00  =  8  j£  2pT-  u 


Corollary  12.  is  not  a  compact  operator  on  C([a,b]). 

Proof.  If  K.  were  compact,  then  Kn,  =Kn-  Knl  would  also  be  compact,  since  Knl  is 

K~  ’  -I  would 


,W 


compact.  Now  G  C([a,b]),  as  V>„(i)  is  bounded  away  from  0.  Hence  ^  _ 

be  compact.  This”  is  impossible,  since  the  identity  operator  I  is  not  compact  on  C(la,bj).  U 
Incidentally,  the  last  proposition  implies  that  if  the  method  of  successive  approximation 
is  applied  to  Equation  (8),  it  will  likely  fail  as  n  increases,  since  PCIL  are  not. uniformly 

We  now  turn  our  attention  to  the  properties  of  Kn  on  C(  * ([<*,&]),  a  €  (Ml-  Unlike 
the  operator  K ,  we  have 

Proposition  13.  C^Qa,  b ])  is  invariant  under  K„  for  a  G  (0, 1]. 

Proof.  It  suffices  to  show  that  gnJ(-) AnJ(-)  G  for  any  cj>  G  ([a,  b]).  If  <f>  G 

C(°,a)([a,fc]),  then  clearly  A n  j<f>  G  b\).  One  can  also  show  that  gn  j  G  C(  '  ’  (K  6])  and 

hence  it  belongs  to  C(°'Q>  ([a,  6])  for  a  G  (0,1].  FinaUy,  <7^  (*)  A„, ,  <£(•)  G  G  (K  b]),  since 
the  latter  is  an  algebra.  U 


To  investigate  the  boundedness  of  Kn  on  C™  ([a,b]),  a  €  (0, 1],  it  suffices  to  consider  the 
individual  components  of  K„ ,  leading  us  to  define  the  following  operators  on  C  ([a,  6]): 

Ln  ■  m  :=  gnJ (i)  A J  =  1,  •  •  •  A- 


Clearly,  FT„  =  £  ^.y  - 

j=i 


9 


„,j  II DO  Halloo)  ^or 


Lemma  14.  L„,  is  bounded  on  CM  ([a,  b}),  a  £  (0, 1],  for  J  =  1,  •  •  • ,  fe.  and  ||I„,  ||„  < 

2 II 9nJ  IL 

Proof.  For  any  <f>  £  C{°'a) {{a,b}),  we  have  \Lni<f>{t)\  =  k,;(0  An,^(0l  <  2U 
all  i  e  [a, 6].  Hence  ||Ln.^|L  <  %JL  M-  For  ^  5  and  *  €  M 

\l^4>(s)  -  =  KAS)  KMS)  ~  ^Al)  KM1) I 

<  I  [s„,jL)  -  5„,i(i)]  A».#^(s)l  +  l^j(<)  [Ar..X5)  _  I 

<  2||^iu  k,»  -  <U*)i  +  Ik, IL  W)  -  <KS)\ 

<  (m\~  k,L  +  lk,IL  \4>L)  I5  -  *1“ 

since  both  gnti  and  <j>  €  C7t#^([a,6|).  Hence  |k,^l«  ^  2M~  k,L  +  l<t>L  ** 

follows  that 

l|£'„j0||„  =  l|b..,^ll-  + 

<  2||sjL  IML  +  2IM-  Is-./ 1-  + 1*- 1-  W- 

<  2|kJIIJI*ll.- 

The  Lemma  is  now  proved.  U 

For  0  <  a  <  1,  one  can  readily  show  that 


2  a°(l  —  a) 

lk,IL  =  ILnjIL  +  k,L  =  y  + 


(l-«) 


Proposition  15.  Kn  is  bounded  on  C{  ’  ’(kb]),  a  6  (0, 1]. 
Proof.  This  follows  directly  from  the  last  Lemma.  U 


While  the  restriction  of  the  operator  if  to  ^“’(k  b] )  allows  it  to  be  defined,  the  restric¬ 
tion  of  the  operator  Kn  to  C'(°,")(k  bl)  does  not  Sain  us  much‘  We  StlU  ^ 


Proposition  16.  K„  is  not  compact  on  &  ’  ([«,&])  f°r  a  e  (°>  11- 

Proof,  The  proof  follows  in  exactly  the  same  manner  as  that  in  the  C([a,b])  case.  U 


Because  K.  are  not  compact  on  CM(M),  we  cannot  make  use  of  the  °f  col¬ 

lectively  compact  operators  to  prove  convergence  of  our  numencal  method.  We  do  have 
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some  type  of  compactness  as  we  will  see  in  the  next  proposition.  However,  this  is  mainly  of 
academic  interest  only. 


Proposition  17.  Kn  :  C 


(o.°) 


([a,  b})  ->  C^0J3)  ([a,  b])  is  compact,  if  0  <  (3  <  a  <  1. 


Proof.  If  we  denote  by  A01'13  the  map  A  :  G(  ’  ’([a,  6])  — >  &  ’([0,6]),  then  K“,f3  —  I  ■ 

Since  Ka-a  is  bounded  and  Ia 3  is  compact,  K*’p  is  compact.  U 


5  Convergence  Theorems  for  Kn 

As  we  cannot  make  use  of  the  theory  of  collectively  compact  operators  to  prove  the  con¬ 
vergence  of  our  numerical  method,  we  resort  to  proving  it  directly.  We  will  prove  some 
pointwise  convergence  properties  of  Kn  after  establishing  several  preliminary  lemmas.  For 
convenience  we  define 


where,  as  before, 


A  <f>{t,s) 


(j,{t)  -  <f>{s) 

g„{t,s)\t- sf 

{5  6  (a,  6]  I  \s-t\<8} 

j  ‘  *  ‘ 

|  0  for  t  =  s 


9sa{p^  5) 


1 

I*-*! 


Lemma  18.  Let  (f>  E  C(°’0)  ([a,  b}),  (3  €  (0,1)  and  a  6  {0,(3].  Then 

|r?(M)l  <  Mh 

for  all  {t,s)  €  [a,h]  x  [a>^]- 

Proof.  This  follows  trivially  from  the  definition  of  C(  '([a, b]). 


Lemma  19.  Let  t  €  [a, 6]  and  <j>  €  C™([a,b}),  (3  €  (0,1].  Then  for  any  e  >  0,  there  exists 
5{e,(j>)  >  0  independent  of  t  such  that 

I  /  5)A^(5>  t)  <^s|  e 
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for  all  8f  <  S. 


Proof.  Let  t  e  [xltxt]  C  [a,b].  Then 


| £  g..(t,s)A^s,t)ds\  <  l’  Mt,s)F*(t,.)\ds 

<11*11,/  lt_  sl  d5  ^  p  (  a  J 


2“*  11*11, 


Thus,  the  required  5  can  be  chosen  as 


r  _  1  (_li- 

2  W-pu 


T  ,  o  ,  ,  c  m  0  a  (o  11.  Then  for  any  e  >  0  there  exist 

Lemma  20.  Let  t  G  [a,  b]  and  <j>  G  6  (Mb  P  M 

AT  (independent  of  t),  0  <  kN<1  <  and  <5  >  0  such  that 

|S„(r):=  £  5„,(r)  A„,</>(r)  1  <  e 

i=fc„,.+1 

for  all  r  G  B(f,  J)  and  for  all  n>N,  where 

a  +  kn,ih„  —  a  +  kN1hN  =:  a, 
a  +  fc„2h„  =  a  +  kN2hN—:b1  and 

£(£,£)  C  KA] 


Proof  Let  e  >  0  be  given.  Let  N  =  «t{»„  %},  where  «,  nre  specified  below. 

For  any  positive  integer  n,  define 

.  _  (  k-.O  ifl^'<fc- 

■  \  k-»»M  = 

.  „  ,.  .  .  .  •  -r  r  r  y-*»  t  G  I  for  a  unique  j  (£),  1  <  J„(0  —  k„- 

Since  [n,h]  is  the  disjoint  nmon  of  {lnJij=i  >  c  Jnj‘io 

4  n  >  L  1 _ [/  £  1  .  It  follows  that  [ai,&i]  2  E*(t)  anc*  ^,i  <  — 

We  define  [aiA]  —  J”U 
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k  for  all  n  >  N.  For  ease  of  presentation,  we  assume  t  E  (aiA)-  (If,  for  example,  t 


t  and  t  ^  a,  we  can  increase  N  by  1,  and  let  [a1}6 1]  - 

* —  1 

following,  we  will  always  assume  r  €  (a,, 6,)  =  (h... l1,. , k„, h„ )  at|l  n  >  N.  Then 

5.W  =  S«W  +  S..,(7-)H-S.,(t) 

where 


III  the 


S„,(r) 

=  E 

^.,(r)  An.,Xr) 

fc„.,  <i<J*  (T) 

5n,(r) 

—  w  ,  q 

n*3n(T) 

.  (r)  A  .  4>{t) 

»v  '  ».»■>>'  ' 

5„,3('r) 

=  E 

9n,AT)  ^J{T) 

We  shall  show  that  each  of  the  above  can  be  made  arbitrarily  small  uniformly  for  r  E  (ai>  h) 
for  n  sufficiently  large.  Indeed, 


provided 


I  S'...  Wl  =  2WC;w)_^(T)l 

<  m\\e\t ...  -Tis 

<  2MX 

<  £ 

~  3 

n>  N2  :=  max  |  [  ^log2  (6  ~  °)P  IMI/9)'l»1j 


We  will  now  derive  a  bound  for  |  5„  ,(t)  |.  If  j*(r)  =  fc,,.  +  1,  then  Snl  (r)  -  0.  Otherwise, 
if  i‘(r)  >  ,  +  1,  then 

|5„,(r)|  <  hn  E 


.  It*  —  T 

<  *.w,  E  E-tj 

kn,l  <J<CM 

=  K  IMI„  E  (X  +  ' 

m=0 


p-\ 
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where  x  :=  r  -  t  .  and  L„  :=  j»  -  -  2.  From  the  definitions  of  and  J», 


we  must  have  ^  <  x  <  3>-  Let  y  :=  Then  |  <  y  <  |  and 

|s„,(r)|  <  <IMU  E  (y  +  ^r1 

m=0 

=  ^ll^Uy^'  +  ^IIC  E  (y  +  m)‘ 

m“  1 

Since  p  e  (0,1]  and  |  <  y  <  |, 

|  S  ,  (r )  |  <  2*~V  Ml  +  K  U\\p  E 

•m  =  l 

From  £(r)  <  fc„.2  and  Ln  =  j»  -  -  2  we  have 

LnK  <  {Kj-KJK 


3-1 


=  ^1  —  fll 


and  so 


|SBtl(r)|  ^  (2'^  +  1)£||«^IU  +  _ai)' 


< 


provided 

where 


n>  Ny  :=  max-fA/i.a,  Ni.b} 


7V1>a  =  max 


{^iog2  (T(fe-a)'W^,1} 

JVlifc  =  max{[^log2  {j^(b-afU\\^j],l} 

Estimating  a  bound  for  |  S„,3(r)  |  is  similar  to  that  for  |  S„,,  (r)  |,  producing  JV3  similar  to 
Nx .  The  Lemma  is  proved  by  choosing  £  =  mm{|i  —  ai|,  \t  M/*  u 

Remark.  It  follows  from  the  proof  of  the  last  lemma  that 

IE  »nJUr) A-»l<6 


j=m1 
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as  long  as  kni  +  1  <m1<m2<  kn2. 

Proposition  21.  For  any  $  €  C'<°'^1  ([«,  &])>  p  G  (0, 1], 

lim  —  KML  ~  0 


P  f  Tet  e  >  0  be  given  By  the  compactness  of  [a,b],  it  suffices  to  show  that  for  any 
Proof,  Let  e  >  U  be  given. .  r  <  e  for  aU  n  >  Nt  and  for  all 

t  in  [a,b],  there  exist  Nt  and  5t  such  that  |(/f  -  Kn)4>\r)\  <  d  for  any 

r  G  B{t,8t).  Now  for  any  7l  and  7,  with  0  <  7j  <  i  -  a,  and  U  S  7,  _ 

ml  and  m2  with  1  <  ^  m2  — 

(K-KU{t)  =  T1(r}a1t-7l)  +  T1(r,t-7l,t  +  72)  +  ^(^f  +  ^’6) 

-  T2(r>»>l,m1)-rj(r,n,m1  +  l,ma)-ra(r>n,ma+  ,  „) 

where 

Tx{t,  x,y)  =  f  5»a(L's)A</)(s,tj  ds 

1  ‘'x 

T2(T,n,m,k)  7  A„1i</’(r) 

j=T71 

» •*-  r :  *•  ■  “T-T™  “  >T  ‘  " 

letting  7l  =  t  -  tNi  ^  tl  >  0  and  -  r«i>Wi ,* 

|T,  (r,  t  —  7l ,  t  +  %)l  —  g 

e 

jT2(r,n,  fcnI  +  1,  fc„,2)l  —  2 

for  aU  n  >  N,  and  r  €  Let  ;7  f  M ]' x  [oT-  7  I’^there  exist  N, 

and  therefore  uniformly  continuous  on  [  ]i  -  J, ,  t  +  i.]  n  K  &1  i  X  l“,  T> ' ' 

and  <$2  >  0  such  that 

|Tl(r,a,t-7l)  -  T2(r,n,  1,^)1 


if/"  l/(r,-)-/(r,C)l*l 

i=i 


e 

<  - 
~  6 
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for  all  n  >  N7  and  r  £  B(t,S2).  Similarly,  there  exist  N3  and  >  0  such  that 

\tAt^  +  72i6)  “  T,iT>n,K,2  +  !A)I 

=  I  E  /  l/(r>«)  -/(r>C,i)]  <**1 

i=fc  ,+i  '»j-i 


< 


for  all  n  >  N3  and  r  £  B(t ,  £,).  The  proposition  is  proved  by  letting  N  =  max{jV,,  AT2,  Nj} 
and  S  =  minj^j,  52,  £,}.  U 

If  follows  immediately  from  the  definition  of  X0  that 


Corollary  22.  For  any  <j>  £  X0, 

lim  ||(tf  -  KM.  =  0 

71— hOO 

If  (f>  e  C{0^]  ([a,  6]),  P  G  (0,  1],  then  we  already  know  that  both  Kef)  and  Kn(f>  G  C{  ’  *  ([«,&]), 
for  ot  G  (0,/?).  One  may  suspect  the  convergence  in  the  last  proposition  is  also  true  in 
CL°'a)  ([a7b}),  i.e.,  lira  \\(K  -  Kn)(f>\\a  =  0  for  0  <  a  <  /?.  Unfortunately,  we  have  not  been 

able  to  prove  it.  However,  for  those  (j)  G  ^([a,  &])  that  satisfy  the  additional  assumption 
in  the  following  proposition,  we  do  have  convergence  in  certain  C{  ’  *([<*,  Z>]). 

Proposition  23.  Let  </>  G  ^([a,  6]),  fi  G  (0,1]  such  that 

\K»m-KJ{s)\<C<\t-sf 

for  some  0  <  C<  P  and  for  some  constant  C ^  independent  of  n,  then 

L'm  UK -KM.  =  0 


for  0  <  a  <  £. 

Proof.  Because  of  the  last  proposition,  it  suffices  to  prove  lim  | (K  —  Kn)(j>\a  =  0.  For 
simplicity  of  notation,  let 

A „K4(t,s)  :=  (K  -  KM1)  -  (*  -  KM>) 
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Then  for  0  <  7  <  7 

\^mt,s)\  =  \Anmt,s)?\KW>s)r 

<  (i km  -  k^)i + 1 Km  -  /cMir  ia„  s)l 

Since  0  €  (^(M),  €  C^(M),  “  «  <  * 

|A„W,5)|  <  {\I<<f>\t+CS \t-s)f\KK<l>{tts)\ 

<  CAt-sfm-KjC 

where  Cx  =  2M(l^l€  +  C*)"-  The  P^Position  f°U°WS  by  *****  7  =  t‘  U‘ 

6  Convergence  Theorems  for  I  7^-n 

.  .  ^  .  A  ■—  T  —  n/K  the  matrix  Bn  defined  below  will 

In  studying  the  properties  of  the  operator  A  -  7  »> 

be  useful.  For  any  7  G  <?««,  6])  with  7(«)  >  0  in  (a> 6)’ 

cnJ(t)  :=  7(*K,sUt)  >  0  on  (a’b)>  j  = 


bn(t)  :=  1  +  £  c„,(t)  >  0  on  la’bl 

i=1 

»u  =  Mc,)<„-UC).  •*«**'• 


Then. 


where 


Bn  is  simply  the  discretized  version  of  A„  in  the  sense  that 

#„(<KC,i));=i . **  =  . 


STS  a.b:"- - 1 . o-  -  vCTify  that  ^ = 1 + A" ' = 

1,  -  •  • ,  K  Thus,  applying  Gerschgorm  Circle  Theorem  [8^  ^  the  eigenvalues  must 

in  the  union  of  the  disks  \z  -  U  <  Ai,  /  ^  1  - 
have  absolute  values  >  1.  Hence,  Bn  is  invertible.  U 
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Proposition  25.  Let  7  £  C([a,b])  and  7 [t)  >  0  on  (a,  b).  Then  for  each  positive  integer  71, 
An  maps  C(J)  1-1  onto  C(J),  where  J  [a,  b\. 


Proof.  An  is  clearly  defined  on  C(J).  For  a  given  x  €  C(J),  we  can  define,  because  of  the 
invertibility  of  Bnl 

X„  :=  (x(ri)>...>x(C1J)t 

Then  it  can  readily  be  shown  that  An(j>  =  where 

««)  =  M‘)"Wi)  +  E^W<l  (is) 

j=l 

Hence  An  is  onto.  If  An<f>  =  0,  then  An<f>{t*  )  =  0,  2  =  l,...,fcn.  Since  is  invertible, 
. )  =  0,  i  —  1, . . .  ,  kn.  Subsequently,  =  0,  since  >  0.  U 

Corollary  26.  yl^1  is  bounded  for  each  n. 

Proof.  Clearly  An  —  I  —  7  is  bounded  on  the  Banach  space  C([a,  t]).  The  boundedness 
of  A  1  is  a  consequence  of  the  Open  Mapping  Theorem  [2].  U 


Other  properties  of  the  matrix  Bn  that  we  will  need  are  contained  in  the  following  lemmas. 


Lemma  27.  Bn  is  irreducible. 

Proof.  Because  Bn  is  a  full  matrix  with  no  non-zero  entries,  it  is  irreducible  [8]. 


Lemma  28.  Bn 1  >  0,  i.e. ,  all  entries  are  positive. 
Proof.  Bn  is  real,  irreducible,  diagonally  dominant  with 


i  ^  j 
i  =  j 


The  Lemma  now  follows  from  a  theorem  in  Varga  [8]  (p.  85).  U 


Lemma  29.  Each  row-sum  of  Bn  is  1. 


(14) 
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Proof.  The  i-th  row- sum 


of  Bn  is  E  bt  j.  From  the  definition  of  biJt  each  is  seen  to  be  one. 


Remarks.  Since  the  entries  bxJ  of  Bn  do  not  all  have  the  same  signs, 


max  |h,,  l  f  !• 


T.ernma  30.  If  each  row-sum  of  a  non-singular  n  X  n  matrix  A  1S  tlien 

has  the  same  property.  «-  ..i  ,, 

Proof.  Let  A,J  =  denote  the  n  X  »  matrix  which  it  identical  to  A  except 


j-th  column  where  it  consists  of  all  ones.  Because 

the  determinant  of  A  is  the  same  as  the  determinant  of  Aj  for  l<j  <  «■  ^  row  SUm 

rf  P  is  the  snm  of  the  cofactors  along  the  i-th  column  of  A  divided  by  the  determinant  of 
A  However,  the  sum  of  the  cofactors  along  the  1-th  column  of  is  jus 
A.  which  is  also  the  determinant  of  A.  Hence  the  i-th  row  sum  of  A  is  one.  U 

CoroUary  31.  \\B?  |L  =  1  for  a11  Positive  integer  n' 

^fT^hi7fonows  from  the  last  lemma  and  the  fact  that  all  entries  in  B  is  positive. 

Proposition  32.  is  uniformly  bounded  on  C(ta>bD' 

Proof.  Let  X  €  C([a,b]).  Then  using  the  notations  in  Equation  (13),  we  have 

Allx{t)  =  K(t)'1\K{t)  +  Gn{t)B^x -1 

Wl"re  C.(*):=(  c„J(t),...,o„.„W),W  =  '=» 

From  the  definition  of  6n(t),  *  follows  im™ediately  that 

\bn(tf\  <  1  and 

<  1 
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for  all  t  G  [a,  &]  and  n  >  0.  Also, 

=  ||X„!L 

<  llxlloo 


It  readily  follows  that 

IKxIL<2||xIL 

for  each  X  in  C([a,b])  and  for  all  n  >  0.  Hence  \\A~J  L  <  2  for  all  ».  U 


Theorem  33,  Let  x  G  <?([«,  6])  and  assume  (I  -  1K)<j>  =  X  has  a  unique  solution  G  X0. 
For  each  positive  integer  n,  let  <j>a  be  the  solution  of 

AJ„  :=  (/  -lKn)<f>n  =  X 


Then 


||<£  -  ^nllco  0 


as  n  — >  oo. 


Proof-  Following  the  standard  arguments,  we  have 

0  '=  {I-7K)<l>-{I-lKn)<l>n 
—  4>  -  <j>n  -  i{K4>  -  Kn4>„) 

=  4>  —  <i>n.  —  -  Kn4>  +  Kn<j>  —  K„<f>n) 

=  t-tn-i&K-KM  +  KM-t*)) 
=  (I  -  7 KM  -  <f>n)  -  l{K  -  Kn)4>  . 

=  An{<(>  -  <!>„)  -  l{K  -  K„)<l> 


Hence 


a:'(K-K.)*. 


II*  -  *„IL  <  cilOl-IK*  - 

Tl,e  theorem  follows  from  the  un.form  boundedness  of  A?  and  the  pointwise  convergence 


of  Kn  to  K.  u 
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7  Conclusion 

In  this  report  we  have  analyzed  the  numerical  solution  of  the  singular  integral  equation 

{I-~tK)<}>  =  x  (15) 

using  the  Nystrom  method  described  in  Section  4.  Here 

K<j)(t)  =  [  1 1-  sr1  [<f>[s)  -  ds 

J  a 

We  studied  the  mapping  properties  of  the  operator  K  and  found  that  the  space  (X0,  ||  ■  |L) 
of  all  uniformly  Hblder  continuous  functions,  despite  not  being  a  Banach  space,  is  a  natural 
setting  to  study  the  unbounded  operator  K,  as  it  (  X0  )  is  invariant  under  K. 

We  also  studied  the  mapping  properties  of  the  numerical  integral  operators  Kn  that  arise 
from  the  Nystrdm  method.  It  is  found  that  Kn  are  bounded  on  C([a,b])  and  (therefore)  on 
X  but  they  are  not  compact  on  C([a,b]).  Nevertheless,  we  proved  a  pointwise  convergence 
theorem  of  K„  to  K  on  (X0,  ||  ■  ||J.  Using  this  and  other  properties  of K  we  proved,  un  er 
appropriate  conditions,  the  convergence  of  the  numerical  solutions  of  the  singular  integral 
equation  (15)  to  its  actual  solution. 
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